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Abstract: Aspects of the molecular-level basis for the function of ethylene glycol 
dimethacrylate and trimethylolproprane trimethacrylate crosslinked methacrylic acid 
copolymers molecularly imprinted with (^-propranolol have been studied using a series 
of all-component and all-atom molecular dynamics studies of the corresponding 
prepolymerization systems. The crosslinking agents were observed to contribute to 
template complexation, and the results were contrasted with previously reported 
template-recognition behavior of the corresponding polymers. Differences in the extent to 
which the two crosslinkers interacted with the functional monomer were identified, and 
correlations were made to polymer-ligand recognition behavior and the results of nuclear 
magnetic resonance spectroscopic studies studies. This study demonstrates the importance 
of considering the functional monomer-crosslinker interaction when designing molecularly 
imprinted polymers, and highlights the often neglected general contribution of crosslinker 
to determining the nature of molecularly imprinted polymer-template selectivity. 
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1. Introduction 

Materials capable of facilitating Angstrom- or nano-scale events such as in chemical catalysis 
or material-biomacromolecular interactions require architectures presenting functionalities and 
molecular-level features that permit high-fidelity molecular discrimination [1]. 

Molecular imprinting [2-8] is a technique for producing highly selective synthetic receptors for 
molecular structures spanning in size from ions to biomacromolecules. The method involves the 
formation of cavities in a synthetic polymer matrix that are of complementary structure and function to 
a template molecule/entity, a ligand for the synthetic receptor. The ability of molecularly imprinted 
polymers (MIPs) to selectively recognize and adsorb the imprinted ligand in the presence of closely 
related chemical species has led to the deployment of these materials in a range of biomedical and 
biotechnological applications, as antibody combining site or receptor binding site mimics. As such, 
MIPs have demonstrated affinities and cross-reactivity profiles comparable to their biological 
counterparts [9-12]. They have been employed as substitutes for biological antibodies in various 
application areas [13], e.g., medical diagnostic, forensic assay development, in solid phase extraction, 
as biosensor recognition elements and in catalysis. Recent examples of the use of MIPs in biological 
contexts [14-16] highlight the potential of molecular imprinting science and technology. The chemical 
and physical stabilities normally associated with molecularly imprinted materials [17] are factors 
further motivating the choice of MIPs as alternatives to biomacromolecules in many areas of application. 

Historically, efforts to understand the molecular-level basis for the molecular imprinting process 
have focused on the interaction of functional monomers and template, as reflected in the development 
of a wide range of novel functional monomers [18-21], a focus that is also reflected in the 
early theoretical studies of MIP systems [22-24]. Throughout the development of the field, there 
has been a general awareness regarding the importance of crosslinking for capturing the functional 
monomer-template complexes and for determining polymer morphology. New crosslinking systems 
have been proposed in a number of studies; e.g., Kempe et al. explored a series of multi-acrylate 
crosslinking agents [25], Spivak et al. [26,27] presented functional monomers incorporating a second 
polymerizable moiety (conversely, crosslinking agents adorned with functionalities often associated 
with functional monomers) and Takeuchi explored the use of disulfide bonds in manipulating 
ligand recognition [28]. In more recent work, Piletsky et al. [29] have examined the possibility of 
orthogonal polymerization strategies where crosslinking of a monomer is accessible by more than one 
approach, Golker et al. [30] have studied the relationships between MlP-template recognition and 
template-crosslinker- functional monomer stoichiometries, and Henschel et al. [31,32] have explored 
the role of crosslinking monomer on MIPs developed for catalysis of the Diels-Alder reaction. Despite 
the acknowledgement of crosslinking as an important feature, little attention has been directed towards 
how the template-crosslinker interaction contributes to the recognition properties of the final material. 
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In a seminal study by Andersson [33], the development of MIPs selective for the P-adrenergic 
antagonist (^-propranolol (SPR) was described. The impact on template recognition of the two 
acrylate-based, crosslinking agents ethylene glycol dimethacrylate (EGDMA) and trimethylolproprane 
trimethacrylate (TRIM) when used in copolymers with the functional monomer methacrylic acid 
(MAA) was examined, as shown in Figure 1. SPR-MIPs have been the subject of a significant number 
of papers on account of interest in new methods for determining the enantiomeric purity of this 
important pharmaceutical and even in the establishment of these materials for use in the detection and 
removal of this compound from drinking water [34-47]. The initial success achieved with SPR-MIPs [33], 
together with propranolol's inherent chirality, a valuable property for the development of reference 
or control studies, as well as its availability in radiolabeled form, have led to it being used as a 
"work-horse" template for fundamental investigations of the molecularly imprinted process, and for 
testing the suitability of MIPs in new application areas. 

Figure 1. Structures of compounds used in this study: (a) the template (^-propranolol 
(SPR); (b) the functional monomer methacrylic acid (MAA); (c) the crosslinking monomer 
ethylene glycol dimethacrylate (EGDMA) and (d) the crosslinking monomer 
trimethylolproprane trimethacrylate (TRIM). 




The desire to improve the efficiency of MIP design and development has motivated a growing 
number of studies aimed at creating a more detailed molecular-level understanding of the molecular 
imprinting process [48-50]. The recent development of full-system, all-atom molecular dynamics 
(MD) studies of MIP prepolymerization systems [51] has provided researchers with the tools for 
examining previously unexplored aspects of the molecular imprinting process [12,30,52-56]. The 
importance of the template-crosslinker interaction for MlP-template recognition has been an important 
observation in some of the systems studied, motivating further study of this aspect of the molecular 
imprinting process. In the present study, we have employed these computational tools to investigate 
the roles of EGDMA and TRIM in the SPR-MIP prepolymerization mixtures, and have drawn 
correlations with template recognition in the corresponding polymers, and with data from ! H-nuclear 
magnetic resonance (NMR) studies of the corresponding prepolymerization mixtures. 

2. Results and Discussion 

2.1. Molecular Dynamics Simulations 

We elected to use full-system all-atom MD-based studies of two polymer systems previously 
described by Andersson [33] (polymers A and C in the original paper) to examine the molecular basis 
for template recognition in the EGDMA- and TRIM-crosslinked MAA copolymers reported in that 
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seminal study. These polymers were prepared using the same solvent, initiator and the same 
template-functional monomer stoichiometry, though with different relative amounts of crosslinker. 
The MD simulations were performed using the same prepolymerization mixture stoichiometries as 
employed in the earlier studies (Table 1). Each system was simulated in quintuplet and interactions 
between species were evaluated using hydrogen bond analyses (Tables 2 and 3). The hydrogen bond 
analyses revealed notable differences between the ensembles of interactions present in the two series 
of prepolymerization mixtures. This is reflected in the relative prevalence of interactions in the 
prepolymerization mixtures (here expressed as percentage of the simulated time) and is summarized 
in Table 2. 

Table 1. Compositions of systems simulated (numbers of molecules present in mixture) a . 



Component 


System A 


System B 


(^-propranolol (SPR) 


10 


10 


Ethylene glycol dimethacrylate (EGDMA) 


398 




Trimethylolproprane trimethacrylate (TRIM) 




80 


Methacrylic acid (MAA) 


80 


80 


Toluene 


1012 


393 


Azobisisobutyronitrile 


6 


2 



a , System stoichiometries are representative of those of the polymers presented by Andersson [33]. 

In the case of the EGDMA copolymers, the crosslinker was observed to be the species present in the 
prepolymerization mixture that contributed most to the complexation of the template (System A, 
Tables 2 and 3). On average, the template is engaged in interactions with the crosslinking monomer for 
roughly 55% of the simulation. This is a phenomenon similar to that observed in MD-studies of other 
EGDMA-MAA copolymers [30,51,52]. The predominance of EGDMA-template interactions was in 
stark contrast to the situation of the functional monomer MAA, which contributes relatively little to the 
template complexation (~9% of time), though has a significant degree of association with EGDMA 
(-62%). Closer analysis revealed that the acidic proton of MAA displays the most stable interactions 
observed (longest average lifetimes, Table 4), in particular with the carbonyl oxygen atoms of 
EGDMA, as well as the highest occupancy values observed in this system (Tables 2 and 3). 

Table 2. Time engaged in hydrogen bond formation (% of simulation) between species in 
Systems A and B. 



System A System B 

Component 

P EGDMA MAA TRIM MAA 

SPR 55.5 8.5 35.8 30.8 

MAA 61.9 n.a. 38.6 n.a. 



n.a., not analyzed. 
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Table 3. Observed average hydrogen bond occupancies a . 









System A 






System B 




Component Atom 


VI A A 

.VI /\/\ 


SPR 


VIA A 
.VI /v/\ 


SPR 






HAA 


H28 


H36 


HAA 


H28 


H36 




\T1 1 

IN 1 1 


2.8 (2.23) b 


Yl.CL. 


n.u. 


12.2 (2.20) 


n.u. 


n.u. 


SPR 

kj 1 IV 




0.1 (0.09) 


II. Li. 


n n 

11. Li. 


0.2 (0.14) 


vi n 

11. Li. 


vi n 




01 s 

V / 1 »y 


Z.O ^ 1 .UJ) J 


11. u. 


vi n 
11. u. 


U.Z yi . /^j 


vi n 


vi n 


MAA 


OAC 


c 

n.u. 


0 1 (0 03") 


0 2 (0 1 4) 


n.u. 


0 2(0 ]]) 


0 1 (0 1 >\) 


had 


n. u. 


1 n (C\ 78A 


1 ^ qa"\ 
i .o (^u.yo j 


n.u. 


A A (\ 9 HA 


7 A (C\ S^A 




04 


0 1 (0 02) 


0.1 (0.07) 


0.1 (0.05) 








EGDMA 


07 
09 
013 


U. 1 v u.uz ^ 

29.9 (2.67) 
31.8(2.02) 


0.1 (0.06) 
11.5 (1.65) 
10.3 (2.05) 


0.0 (0.02) 
16.4 (4.11) 
17.0 (3.20) 










015 








13.5 (1.15) 


3.8(1.40) 


8.2 (3.39) 




022 








12.3 (1.68) 


3.8(2.10) 


8.1 (5.77) 


TRIM 


023 








12.7(1.82) 


3.6(1.38) 


8.4 (3.41) 


04 








0.0 (0.00) 


0.0 (0.00) 


0.0 (0.00) 



O10 - - - 0.0(0.00) 0.0(0.00) 0.0(0.00) 
017 - - - 0.0 (0.00) 0.0 (0.00) 0.0 (0.00) 

a , The values were calculated by summation of all observed occupancies (in percentage of simulation time) 
and for each analysed interaction and division of this sum by the number of reference molecules for each 
system. For all interactions involving the template (SPR) the occupancy values are calculated as "per 
template". For all monomer-monomer contacts, the values were averaged against the number of functional 
monomers (MAA) being present at constant amounts in the evaluated systems. These values where then, 
again, summarized from each of the quintuplet simulations and a total average calculated from these results; 

b , Values in brackets are the standard deviations for the average occupancies from quintuplet simulations; 

c , n.a., not analyzed. 

When examining the TRIM-crosslinked polymers (System B), it was found that the extent of 
interaction of the crosslinker with the template was lower than in the EGDMA-system (Tables 2 and 3). 
The extent of interaction of TRIM with the functional monomer was notably lower and in clear 
contrast to the situation in the EGDMA system. Further, the SPR-MAA interaction was considerably 
more prevalent in the TRIM system than in the case with EGDMA (-31% compared with -9%, 
respectively). We also noted that from the perspective of the template, the total extent of interaction 
with monomers (crosslinker or MAA) over the 10 ns were quite similar, 64.0% and 66.6% for Systems A 
and B, respectively. 

Initially, we attributed these differences in the extents of interactions to differences in the 
SPR: crosslinker and MAA: crosslinker stoichiometries employed in the original study, i.e., 
40 EGDMA molecules per template (System A) and 8 TRIM molecules per template (System B), 
corresponding to 80 carbonyl oxygens per template in System A and 24 carbonyl oxygen atoms per 
template in System B. Upon closer examination, we considered that the relative differences observed 
in the MD data (Tables 2 and 3) could be dependent upon the number of crosslinker-derived hydrogen 
bond accepting atoms potentially available for interaction with the template, however the ~3-fold 
lower total number of potential H-bond accepting carbonyl oxygens in the TRIM-system cannot alone 
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account for the ~0. 3-fold lower degree of interaction of the template with this crosslinker, as compared 
to EGDMA. Furthermore, in the case of the crosslinker-MAA interactions, the extent of MAA's 
interaction of MAA with the respective crosslinkers did not directly correlate with the differences in 
numbers of hydrogen bond acceptors available. Collectively, this suggested that despite the similarities 
in the functionalities present in these crosslinkers, there might be differences of a more fundamental 
character between these two cross-linking monomers when used in these two polymer systems. 



Table 4. Observed average hydrogen bond lifetimes a . 









System A 






System C 




Component 


Atom 


MAA 


SPR 


MAA 


SPR 






HAA 


H28 


H36 


HAA 


H28 


H36 




Nil 


1.7 (0.30) b 


n.a. 


n.a. 


1.9 (0.04) 


n.a. 


n.a. 


SPR 


06 


0.8(0.13) 


n.a. 


n.a. 


0.7 (0.03) 


n.a. 


n.a. 




015 


1.7 (0.13) 


n.a. 


n.a. 


1.7 (0.14) 


n.a. 


n.a. 


MAA 


OAC 


n.a. c 


0.6 (0.03) 


0.7(0.10) 


n.a. 


0.6 (0.02) 


0.8 (0.06) 


OAD 


n.a. 


0.8 (0.05) 


1.4 (0.47) 


n.a. 


0.8 (0.03) 


1.6 (0.11) 




04 


0.7 (0.03) 


0.6 (0.06) 


0.6 (0.03) 








EGDMA 


07 
09 


0.7 (0.03) 
3.4 (0.09) 


0.6 (0.05) 
0.8 (0.05) 


0.6 (0.05) 
1.5 (0.15) 










013 


3.3 (0.13) 


0.8 (0.04) 


1.7 (0.10) 










015 








3.3 (0.06) 


0.9 (0.06) 


1.7 (0.26) 




022 








3.1 (0.20) 


0.8 (0.03) 


1.6 (0.26) 


TRIM 


023 








3.1 (0.15) 


0.8 (0.04) 


1.7 (0.26) 


04 
O10 
017 








0.6 (0.03) 
0.6 (0.03) 
0.5 (0.03) 


0.0 (0.00) 
0.0 (0.00) 
0.0 (0.00) 


0.0 (0.00) 
0.0 (0.00) 
0.0 (0.00) 


a , All time 


values 


are presented 


in picoseconds. 


The avera^ 


*e interaction 


lifetimes were 


calculated 



by summation of all hydrogen bond event lifetimes and divided with the observed number of events in 
each quintuplet system. These average lifetimes were then summarized and averaged as the occupancies 
in Table 3. b , Values in brackets are standard deviations of the average values presented in this table. 
c , n.a., not analyzed. 

2.2. Correlations with Polymer-Template Binding and NMR Studies 

To examine this further, comparisons with polymer-template recognition data were undertaken. 
In the original description of the binding of template by the EGDMA- and TRIM-based imprinted 
polymers [33], the IC 5 o values for binding in toluene were determined by radioligand binding 
studies to be 0.29 and 0.16 uM, respectively. Interestingly, the difference between these systems is 
reflected in the difference in the extent of MAA-SPR interactions observed in the MD studies of the 
prepolymerization mixtures of these polymers (stronger binding/lower IC 5 o and greater degree of 
MAA-template interaction). This suggested that the relatively limited MAA-TRIM interaction 
(in relation to the number of hydrogen-bond acceptors in each system) could provide a basis for the 
relative prevalence of the SPR-MAA interaction in System B. The basis for the lower extent of 
MAA-TRIM interaction per H-bond acceptor, as compared to the EGDMA system, was attributed to 
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the comparatively greater steric crowding of hydrogen-bond acceptors of TRIM, relative to the 
situation in the EGDMA-based system. In other words, the third carbonyl of TRIM does not offer the 
same potential for interaction with the functional monomer or template as the first two, or those of 
EGDMA. This situation, in turn, leaves a greater proportion of the MAA available for relatively 
stronger interactions with the template, e.g., with the amine. While the degree to which the 
crosslinking afforded by the crosslinking monomers contributes to polymer-ligand recognition is not 
addressed here, it is the subject of ongoing theoretical and experimental studies in our laboratory. To 
the best of our knowledge, this is the first study in which the interplay between crosslinker and 
functional monomer can be seen to influence the degree of functional monomer-template interaction. 

To provide further insight, a series of ! H-NMR studies was undertaken where mixtures of SPR and 
MAA (the template and functional monomer had the same relative molar stoichiometrics in both 
systems) were titrated with the crosslinkers, as shown in Figure 2. Using the methyl resonance of the 
isopropyl group of the template as a diagnostic tool and a one-site binding model, apparent K& values 
of 696 ± 57 mM (standard error of the mean) and 641 ± 103 mM were determined for the 
crosslinker-SPR/MAA interaction in the EGDMA and TRIM systems, respectively. That the chemical 
shift changes of the template's two equivalent methyl groups demonstrated no significant difference 
between the crosslinkers provides support for the interplay mechanism suggested by the MD-studies, 
whereby the competition of crosslinker and template for interaction with the functional monomer in 
both systems provides similar total degrees of complexation. Furthermore, the NMR data demonstrate 
the capacities of the two crosslinking agents to compete, albeit weakly, for coordination of the 
template in the presence of functional monomer. How the impact of crosslinker on the availability of 
functional monomer influences non-specific binding, and even how the degree of crosslinking may 
influence recognition are topics currently being addressed in our laboratory. This study also highlights 
the potential of functionalized cross-linking agents, such as those developed by Spivak et al. [26-27], 
for creating more homogeneous polymer recognition sites. 

Figure 2. ^-NMR titration plots where constant concentrations of propranolol and 
methacrylic acid were titrated with increasing amounts of crosslinking monomer, in (A) 
TRIM and in (B) EGDMA. The isopropyl methyl resonances of (^-propranolol were 
studied (Figures 1 and 3). 
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3. Experimental Section 

3.1. Chemicals 

(^-propranolol hydrochloride, EGDMA, TRIM and (ig-toluene were obtained from Sigma-Aldrich 
(Munich, Germany). Toluene was purchased from Merck (Solna, Sweden). MAA was distilled under 
reduced pressure and kept at -20 °C until used, (^-propranolol was extracted from (^-propranolol 
hydrochloride in a basic aqueous solution (pH ~ 8) using diethyl ether as organic phase. EGDMA and 
TRIM were each extracted three times with a mixture of 75 mL 0.1 M NaOH and 25 mL saturated 
NaCl followed by 25 mL of saturated NaCl. The washed substances were dried over MgS0 4 and 
passed through AIO3 prior to use. 

3.2. Molecular Dynamics Simulations 

All-atom MD simulations were conducted using the AMBER (v. 10.0 UCSF, San Francisco, CA, 
USA) suite of programs [57,58] using a strategy previously described by Karlsson et al. [51] and 
developed in subsequent studies [30,52,55]. Implemented force field(s) in these simulations were the 
AMBER (FF03) [59] and the general amber force field (GAFF) [60]. Simulated systems were initially 
constructed using the PACKMOL software [61,62] to obtain random starting geometries, as presented 
by Olsson et al. [52] The compositions of the simulated all-component prepolymerization mixtures, as 
well as equilibration and production run data, are summarized in Table 1 . All systems were simulated 
(in quintuplet) for 10 ns of recorded trajectory data. 

Production run trajectory datasets were analyzed using HBOND routine available in the PTRAJ 
module included in AmberTools (v. 1.3, UCSF, San Francisco, CA, USA) [57]. All hydrogen bond 
(HBOND) interactions were extracted from the trajectories using a distance and angle cut-off of 3.0 A 
and 60°, respectively. The structures of the analyzed molecular species and the interacting atoms 
potentially participating in hydrogen bond interactions are presented in Figure 3. 

Figure 3. Structures of compounds used in this study with molecular abbreviations and 
atomic labels implemented in performed calculations and evaluations of modeled systems. 




Metliaciylic acid (022)0 0(023) 

(MAA) Tiimethylolproprane trimethacrylate 

(TRIM) 
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3.3. 'H-NMR Studies 

All NMR studies were performed using dg-toluene as solvent on a Varian 500 MHz instrument 
(PaloAlto, CA, USA). Spectra were acquired at 293 ± 1 K using the Varian L700 Pulsed Gradient 
Driver steered with VNMR 6. IB software as supplied by the manufacturer. 

Titration studies followed a previously described general protocol [63]: 17.7 mM (^-propranolol 
solutions containing 142 mM MAA in dg-toluene were titrated using a 2 M EGDMA or TRIM 
solutions containing 17.7 mM propranolol and 142 mM MAA. Between 20 and 32 measurement points 
were used for each study. 

4. Conclusions 

Aspects of the molecular-level basis for the function of the crosslinking monomers EGDMA and 
TRIM in MAA copolymers imprinted with (^-propranolol were studied using a series of 
all-component and all-atom MD simulations of the corresponding prepolymerization systems. The 
crosslinking agents were observed to make significant contributions to template complexation, though 
the MD studies reveal considerable differences in the balance of the ensembles of equilibria present in 
the polymerization systems. Our results indicate that differences in the extent that the crosslinkers 
interact with the functional monomer MAA contribute to its availability for interaction with the 
template, an effect not previously discussed in the literature. This study highlights the often neglected 
contribution of crosslinker in determining the nature of MlP-template selectivity. The results have 
correlated with the results of previously reported template-recognition studies of the corresponding 
polymers and with NMR studies of prepolymerization mixtures provided support for the conclusions 
drawn from the MD-studies. Finally, this study highlights the utility of all-atom full-system molecular 
dynamics studies as a diagnostic tool for the evaluation of molecular imprinting systems. 

Acknowledgments 

The authors acknowledge the financial support provided by the Swedish Knowledge Foundation 
(KKS, grant 2010-0223), Swedish Research Council (VR, grant 2006-6041) and Linnaeus University. 

Author Contributions 

Molecular dynamics studies and data analyses were performed by Siamak Shoravi, 
Gustaf D. Olsson and Bjorn C. G. Karlsson. NMR studies were performed by Siamak Shoravi and 
Gustaf D. Olsson. All authors contributed to the drafting of the manuscript. 

Conflicts of Interest 

The authors declare no conflict of interest. 
References 

1. Zhang, S.G. Fabrication of novel biomaterials through molecular self-assembly. Nat. Biotechnol. 
2003,27, 1171-1178. 



Int. J. Mol. Sci. 2014, 15 



10631 



2. Sellergren, B. Molecularly Imprinted Polymers: Man-Made Mimics of Antibodies and Their 
Applications in Analytical Chemistry; Elsevier: Amsterdam, The Netherlands, 2001. 

3. Alexander, C; Andersson, H.S.; Andersson, L.L; Ansell, R.J.; Kirsch, N.; Nicholls, LA.; 
O'Mahony, J.; Whitcombe, M.J. Molecular imprinting science and technology: A survey of the 
literature for the years up to and including 2003. J. Mol. Recognit. 2006, 19, 106-180. 

4. Whitcombe, M.J.; Kirsch, N.; Nicholls, LA. Molecular imprinting science and technology: 
A survey of the literature for the years 2004-2011. J. Mol. Recognit. 2014, 2014, 
doi:10.1002/jmr.2347. 

5. Wulff, G. Molecular imprinting in cross-linked materials with the aid of molecular templates — A 
way towards artificial antibodies. Angew. Chem. Int. Ed. Engl. 1995, 34, 1812-1832. 

6. Haupt, K.; Mosbach, K. Molecularly imprinted polymers and their use in biomimetic sensors. 
Chem. Rev. 2000, 100, 2495-2504. 

7. Whitcombe, M.J.; Chianella, L; Larcombe, L.; Piletsky, S.A.; Noble, J.; Porter, R.; Horgan, A. 
The rational development of molecularly imprinted polymer-based sensors for protein detection. 
Chem. Soc. Rev. 2011, 40, 1547-1571. 

8. Komiyama, M.; Takeuchi, T.; Mukawa, T.; Asanuma, H. Molecular Imprinting; Wiley- VCH: 
Weinheim, Germany, 2004. 

9. Vlatakis, G.; Andersson, L.L; Miiller, R.; Mosbach, K. Drug assay using antibody mimics made 
by molecular imprinting. Nature 1993, 361, 645-647. 

10. Andersson, L.L; Miiller, R.; Vlatakis, G.; Mosbach, K. Mimics of the binding sites of opioid 
receptors obtained by molecular imprinting of enkephalin and morphine. Proc. Natl. Acad. 
Sci. USA. 1995, 92, 4788-4792. 

11. Berglund, L; Nicholls, LA.; Lindbladh, C; Mosbach, K. Recognition in molecularly imprinted 
polymer a2-adrenoreceptor mimics. Bioorg. Med. Chem. Lett. 1996, 6, 2237-2242. 

12. Schillinger, E.; Moder, M.; Olsson, G.D.; Nicholls, LA.; Sellergren, B. An artificial estrogen 
receptor through combinatorial imprinting. Chem. Eur. J. 2012, 18, 14773-14783. 

13. Sellergren, B.; Allender, C.J. Molecularly imprinted polymers: A bridge to advanced drug 
delivery. Adv. DrugDeliv. Rev. 2005, 57, 1733-1741. 

14. Hoshino, Y.; Koide, H.; Urakami, T.; Kanazawa, H.; Kodama, T.; Oku, N.; Shea, K.J. Recognition, 
neutralization, and clearance of target peptides in the bloodstream of living mice by molecularly 
imprinted polymer nanoparticles: A plastic antibody. J. Am. Chem. Soc. 2010, 132, 6644—6645. 

15. Hoshino, Y.; Koide, H.; Furuya, K.; Haberaecker, W.W.; Lee, S.-H.; Kodama, T.; Kanazawa, H.; 
Oku, N.; Shea, K.J. The rational design of a synthetic polymer nanoparticle that neutralizes a toxic 
peptide in vivo. Proc. Natl. Acad. Sci. USA 2012, 109, 33-38. 

16. Piletska, E.V.; Stavroulakis, G.; Larcombe, L.D.; Whitcombe, M.J.; Sharma, A.; Primrose, S.; 
Robinson, G.K.; Piletsky, S.A. Passive control of quorum sensing: Prevention of pseudomonas 
aeruginosa biofilm formation by imprinted polymers. Biomacromolecules 2011, 12, 1067-1071. 

17. Svenson, J.; Nicholls, LA. On the thermal and chemical stability of molecularly imprinted 
polymers. Anal. Chim. Acta 2001, 435, 19-24. 

18. Tanabe, K.; Takeuchi, T.; Matsui, J.; Ikebukuro, K.; Yano, K.; Karube, I. Recognition of 
barbiturates in molecularly imprinted copolymers using multiple hydrogen-bonding. J. Chem. 
Soc. Chem. Commun. 1995, 1995, 2303-2304. 



Int. J. Mol. Sci. 2014, 15 



10632 



19. Piletsky, S.A.; Andersson, H.S.; Nicholls, LA. Combined hydrophobic and electrostatic 
interaction-based recognition in molecularly imprinted polymers. Macromolecules 1999, 32, 633-636. 

20. Wulff, G.; Schonfeld, R. Polymerizable amidines — Adhesion mediators and binding sites for 
molecular imprinting. Adv. Mater. 1998, 10, 957-959. 

21. Hall, A.J.; Manesiotis, P.; Emgenbroich, M.; Quaglia, M.; de Lorenzi, E.; Sellergren, B. 
Urea host monomers for stoichiometric molecular imprinting of oxyanions. J. Org. Chem. 2005, 
70, 1732-1736. 

22. Subrahmanyam, S.; Piletsky, S.A.; Piletska, E.V.; Chen, B.; Karim, K.; Turner, A.P. 
"Bite-and-Switch" approach using computationally designed molecularly imprinted polymers for 
sensing of creatinine. Biosens. Bioelectron. 2001, 16, 631-637. 

23. Chianella, I.; Lotierzo, M.; Piletsky, S.A.; Tothill, I.E.; Chen, B.; Karim, K.; Turner, A.P.F. 
Rational design of a polymer specific for microcystin-LR using a computational approach. 
Anal. Chem. 2002, 74, 1288-1293. 

24. Chianella, I.; Karim, K.; Piletska, E.V.; Preston, C; Piletsky, S.A. Computational design and 
synthesis of molecularly imprinted polymers with high binding capacity for pharmaceutical 
applications-model case: Adsorbent for abacavir. Anal. Chim. Acta 2006, 559, 73-78. 

25. Kempe, M. Antibody-mimicking polymers as chiral stationary phases in HPLC. Anal. Chem. 
1996, 68, 1948-1953. 

26. Sibrian- Vazquez, M.; Spivak, D.A. Enhanced enantioselectivity of molecularly imprinted 
polymers formulated with novel cross-linking monomers. Macromolecules 2003, 36, 5105-51 13. 

27. Sibrian- Vazquez, M.; Spivak, D.A. Molecular imprinting made easy. J. Am. Chem. Soc. 2004, 
126, 7827-7833. 

28. Takeda, K.; Kuwahara, A.; Ohmori, K.; Takeuchi, T. Molecularly imprinted tunable binding sites 
based on conjugated prosthetic groups and ion-paired cofactors. J. Am. Chem. Soc. 2009, 131, 
8833-8838. 

29. Lakshmi, D.; Whitcombe, M.J.; Davis, F.; Chianella, I.; Piletska, E.V.; Guerreiro, A.; 
Subrahmanyam, S.; Brito, P.S.; Fowler, S.A.; Piletsky, S.A. Chimeric polymers formed from a 
monomer capable of free radical, oxidative and electrochemical polymerisation. Chem. Commun. 
2009, 2009, 2759-2761. 

30. Golker, K.; Karlsson, B.C.G.; Olsson, G.D.; Rosengren, A.M.; Nicholls, LA. Influence of 
composition and morphology on template recognition in molecularly imprinted polymers. 
Macromolecules 2013, 46, 1408-1414. 

31. Kirsch, N.; Hedin-Dahlstrom, J.; Henschel, FL; Whitcombe, M.J.; Wikman, S.; Nicholls, LA. 
Molecularly imprinted polymer catalysis of a Diels-Alder reaction. J. Mol. Catal. B Enzym. 2009, 
58, 110-117. 

32. Henschel, H.; Kirsch, N.; Hedin-Dahlstrom, J.; Whitcombe, M.J.; Wikman, S.; Nicholls, LA. 
Effect of the cross-linker on the general performance and temperature dependent behaviour of a 
molecularly imprinted polymer catalyst of a Diels-Alder reaction. J. Mol. Catal. B Enzym. 2011, 
72, 199-205. 

33. Andersson, L.I. Application of molecular imprinting to the development of aqueous buffer and 
organic solvent based radioligand binding assays for (^-propranolol. Anal. Chem. 1996, 68, 1 1 1-1 17. 



Int. J. Mol. Sci. 2014, 15 



10633 



34. Nilsson, S.; Schweitz, L.; Petersson, M. Three approaches to enantiomer separation of 
P-adrenergic antagonists by capillary electrochromatography. Electrophoresis 1997, 18, 884-890. 

35. Haginaka, J.; Sakai, Y.; Narimatsu, S. Uniform-sized molecularly imprinted polymer material for 
propranolol recognition of propranolol and its metabolites. Anal. Sci. 1998, 14, 823-826. 

36. Mullett, W.M.; Martin, P.; Pawliszyn, J. In- tube molecularly imprinted polymer solid-phase 
microextraction for the selective determination of propranolol. Anal. Chem. 2001, 73, 2383-2389. 

37. Schweitz, L.; Andersson, L.L; Nilsson, S. Rapid electrochromatographic enantiomer separations 
on short molecularly imprinted polymer monoliths. Anal. Chim. Acta 2001, 435, 

38. Fireman-Shoresh, S.; Popov, I.; Avnir, D.; Marx, S. Enantioselective, chirally templated sol-gel 
thin films. J. Am. Chem. Soc. 2005, 127, 2650-2655. 

39. Schmidt, R.H.; Haupt, K. Molecularly imprinted polymer films with binding properties enhanced 
by the reaction-induced phase separation of a sacrificial polymeric porogen. Chem. Mater. 2005, 
17, 1007-1016. 

40. Schmidt, R.H.; Belmont, A.S.; Haupt, K. Porogen formulations for obtaining molecularly 
imprinted polymers with optimized binding properties. Anal. Chim. Acta 2005, 542, 1 18-124. 

41. Castell, O.K.; Allender, C.J.; Barrow, D.A. Novel biphasic separations utilising highly selective 
molecularly imprinted polymers as biorecognition solvent extraction agents. Biosens. Bioelectron. 
2006, 22, 526-533. 

42. Perez-Moral, N.; Mayes, A.G. Direct rapid synthesis of MIP beads in SPE cartridges. 
Biosens. Bioelectron. 2006, 21, 1798-1803. 

43. Ye, L.; Yoshimatsu, K.; Kolodziej, D.; da Cruz Francisco, J.; Dey, E.S. Preparation of 
molecularly imprinted polymers in supercritical carbon dioxide. J. Appl. Polym. Sci. 2006, 102, 
2863-2867. 

44. Reimhult, K.; Yoshimatsu, K.; Risveden, K.; Chen, S.; Ye, L.; Krozer, A. Characterization of 
QCM sensor surfaces coated with molecularly imprinted nanoparticles. Biosens. Bioelectron. 

2008, 23, 1908-1914. 

45. Bompart, M.; Gheber, L.A.; de Wilde, Y.; Haupt, K. Direct detection of analyte binding to single 
molecularly imprinted polymer particles by confocal Raman spectroscopy. Biosens. Bioelectron. 

2009, 25, 568-571. 

46. Long, Y.; Philip, J.Y.N. ; Schillen, K.; Liu, F.; Ye, L. Insight into molecular imprinting 
in precipitation polymerization systems using solution NMR and dynamic light scattering. 
J. Mol. Recognit. 2011, 24, 619-630. 

47. Xu, C; Ye, L. Clickable molecularly imprinted nanoparticles. Chem. Commun. 2011, 47, 6096-6098. 

48. Nicholls, I.A.; Andersson, H.S.; Charlton, C; Henschel, H.; Karlsson, B.C.G.; Karlsson, J.G.; 
O'Mahony, J.; Rosengren, A.M.; Rosengren, J.K.; Wikman, S. Theoretical and computational 
strategies for rational molecularly imprinted polymer design. Biosens. Bioelectron. 2009, 25, 543-552. 

49. Nicholls, I.A.; Andersson, H.S.; Golker, K.; Henschel, H.; Karlsson, B.C.G.; Olsson, G.D.; 
Rosengren, A.M.; Shoravi, S.; Suriyanarayanan, S.; Wiklander, J.G.; et al. Rational design of 
biomimetic molecularly imprinted materials: Theoretical and computational strategies for guiding 
nanoscale structured polymer development. Anal. Bioanal. Chem. 2011, 400, 1771-1786. 

50. Nicholls, I.A.; Karlsson, B.C.G.; Olsson, G.D.; Rosengren, A.M. Computational strategies for the 
design and study of molecularly imprinted materials. Ind. Eng. Chem. Res. 2013, 52, 13900-13909. 



Int. J. Mol. Sci. 2014, 15 



10634 



51. Karlsson, B.C.G.; O'Mahony, J.; Karlsson, J.G.; Bengtsson, H.; Eriksson, L.A.; Nicholls, LA. 
Structure and dynamics of monomer-template complexation: An explanation for molecularly 
imprinted polymer recognition site heterogeneity. J. Am. Chem. Soc. 2009, 131, 13297-13304. 

52. Olsson, G.D.; Karlsson, B.C.G.; Shoravi, S.; Wiklander, J.G.; Nicholls, LA. Mechanisms 
underlying molecularly imprinted polymer molecular memory and the role of crosslinker: 
Resolving debate on the nature of template recognition in phenylalanine anilide imprinted 
polymers. J. Mol. Recognit. 2012, 25, 69-73. 

53. Atta, N.F.; Hamed, M.M.; Abdel-Mageed, A.M. Computational investigation and synthesis of 
a sol-gel imprinted material for sensing application of some biologically active molecules. 
Anal. Chim. Acta 2010, 667, 63-70. 

54. Karlsson, B.C.G.; Rosengren, A.M.; Naslund, L; Andersson, P.O.; Nicholls, LA. Synthetic human 
serum albumin sudlow I binding site mimics. J. Med. Chem. 2010, 53, 7932-7937. 

55. Olsson, G.D.; Karlsson, B.C.G.; Schillinger, E.; Sellergren, B.; Nicholls, LA. Theoretical studies 
of 17-P-estradiol-imprinted prepolymerization mixtures: Insights concerning the roles of 
cross-linking and functional monomers in template complexation and polymerization. Ind. Eng. 
Chem. Res. 2013, 52, 13965-13970. 

56. O'Mahony, J.; Karlsson, B.C.G.; Mizaikoff, B.; Nicholls, LA. Correlated theoretical, 
spectroscopic and X-ray crystallographic studies of a non-covalent molecularly imprinted 
polymerisation system. Analyst 2007, 132, 1 161-1 168. 

57. Case, D.A.; Cheatham, T.E.; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M.; Onufriev, A.; 
Simmerling, C; Wang, B.; Woods, R.J. The Amber biomolecular simulation programs. 
J. Comput. Chem. 2005, 26, 1668-1688. 

58. Case, D.A.; Darden, T.A.; Cheatham, T.E.I.; Simmerling, C.L.; Wang, J.; Duke, R.E.; Luo, R.; 
Crowley, M.; Walker, R.C.; Zhang, W.; et al. Amber 10; University of California: San Francisco, 
CA, USA, 2008. 

59. Duan, Y.; Wu, C; Chowdhury, S.; Lee, M.C.; Xiong, G.; Zhang, W.; Yang, R.; Cieplak, P.; Luo, R.; 
Lee, T.; et al. A point-charge force field for molecular mechanics simulations of proteins based on 
condensed-phase quantum mechanical calculations. J. Comput. Chem. 2003,24, 1999-2012. 

60. Wang, L; Wolf, R.M.; Caldwell, J.W.; Kollman, P. A.; Case, D.A. Development and testing of a 
general amber force field. J. Comput. Chem. 2004, 25, 1 157-1 174. 

61. Martinez, J.M.; Martinez, L. Packing optimization for automated generation of complex 
system's initial configurations for molecular dynamics and docking. J. Comput. Chem. 2003, 24, 
819-825. 

62. Martinez, L.; Andrade, R.; Birgin, E.G.; Martinez, J.M. PACKMOL: A package for building 
initial configurations for molecular dynamics simulations. J. Comput. Chem. 2009, 30, 2157-2164. 

63. Svenson, L; Andersson, H.S.; Piletsky, S.A.; Nicholls, LA. Spectroscopic studies of the molecular 
imprinting self-assembly process. J. Mol. Recognit. 1998, 11, 83-86. 

© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http ://creativecommons. org/licenses/by/3 . 0/) . 



